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Abstract 

We investigate the minimal supergravity (mSUGRA) signals at the LHC in the context of 
supercritical string cosmology (SSC). In this theory, the presence of a time dependent dilaton 
provides us with a smoothly evolving dark energy and modifies the dark matter allowed region 
of the mSUGRA model with standard cosmology. Such a dilaton dilutes the supersymmetric 
dark matter density (of neutralinos) by a factor 0(10) and consequently the regions with too 
much dark matter in the standard scenario are allowed in the SSC. The final states expected 
at the LHC in this scenario, unlike the standard scenario, consist of Z bosons, Higgs bosons, 
and/or high energy taus. We show how to characterize these final states and determine the 
model parameters. Using these parameters, we determine the dark matter content and the 
neutralino-proton cross section. All these techniques can also be applied to determine model 
parameters in SSC models with different SUSY breaking scenario s. 



1 Introduction 



The recent WMAP data [T] has determined the content of the universe very precisely. The 
dark matter and dark energy compose about 23% and 73% of the total energy density of the 
universe, respectively. 

The origin of dark matter can be explained in supersymmetry (SUSY) models, where the 

lightest SUSY particle, the neutralino (in most SUSY models) |2], is the dark matter candidate. 

SUSY combined with supergravity grand unification (SUGRA GUT) [3], resolves a number of 

the problems inherent in the standard model (SM). The SUGRA GUT model not only solves 

the gauge hierarchy problem and predicts grand unification at the GUT scale Mq ~ 10 16 GeV 

but also allows for the spontaneous breaking of SUGRA at the Mq scale in a hidden sector, 

leading to an array of soft breaking masses. The renormalization group equations (RGEs) 

then show that this breaking of SUGRA leads naturally to the breaking of SU(2)xU(l) of the 

SM at the electroweak scale [I]. SUSY breaking masses around a TeV for most of the SUSY 

parameter space are allowed by other experimental constraints. It is also very interesting to note 

that achieving the WMAP relic density requires the annihilation cross section of the lightest 

neutralino in these SUSY models to be of order 1 pb with M~o ~ (9(100 GeV). Such a mass 

X.1 

scale is reachable at the LHC. 

The origin of dark energy is not well understood. The simplest proposal is to add a cos- 
mological constant in Einstein's equation. However, the reason why the dark matter content is 
comparable to the dark energy content at the present time remains a puzzle. Another proposal 
is that a quintessence scalar field is responsible for dark energy [5|. However, this requires the 
field to have a very small mass and is not well motivated in particle physics models. In the 
context of string theory, the dilaton can play the role of dark energy [61 [7] . One also finds pro- 
posals which involve, for example, modifications to General Relativity, Braneworld scenarios, 
or Topological defects, which are invoked to explain this fundamental issue. 

In this paper, we will investigate experimental signatures of SUSY as consequences of a 
rolling dilaton in the Q-cosmology scenario [6] which offers an alternative framework that 
establishes the Supercritical (or non-critical) String Cosmology (or SSC). In the SSC framework, 
the dark energy has two components: One component arises from the dilaton, 0, and the other 
arises from the Q 2 which is associated with the central charge deficit. Both Q and the dilaton 
have time dependent pieces. It was shown that the SSC scenario [8] is consistent with the 
smoothly evolving dark energy at least for the last ten billion years (0 < z < 1.6), in accordance 
with the very recent observations on supernovae [9]. 

The presence of this time dependent dilaton affects the relic density calculation since it 
modifies the Boltzman equation in the following way: % + 3Hn + (av)(n 2 — n e q 2 ) — (fin = 
The relic density is given by 

nh 2 = Rx (nh 2 ) (i) 

where R ~ exp[ ((fiH' 1 /x)dx] and (flh 2 )o denotes the relic density that is obtained by 
ordinary cosmology. It is possible to determine R by solving for <fi from the field equations for 
this SSC scenario. The value of R is about 0.1 in order to satisfy the recent observation of the 
evolution of dark energy in the range < z < 1.6. This new factor changes the profile of dark 
matter allowed region in SUSY models. 
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To investigate the SUSY signatures we use the minimal SUGRA (mSUGRA) model and 
calculate the dark matter content in the context of the SSC framework. We note that the 
low energy limit of string theory is certainly much more complicated than mSUGRA, and 
there are many different effective theories, depending on the details of compactification and 
SUSY breaking [10]. The relevant dark matter phenomenological analyses are highly model 
dependent [llj. In some cases, such as the orbifold-compactified heterotic models [10], there 
might be situations in which the couplings of matter with stabilized dilatons lead to non-thermal 
dark matter, thus leading to completely different phenomenology. 

However, the SSC framework is characterized by a non-stabilized dilaton which runs in 
cosmic time [8]. In this context, it is possible to have thermalization of weakly interacting dark 
matter, such as the mSUGRA lightest neutralino (xi) which couples to the dilaton. In this 
sense, the mSUGRA framework provides a sufficiently non-trivial and generic pilot study of 
the novel effects the running dilaton has on the abundance of thermal dark matter relics. 

The mSUGRA parameters are the universal scalar mass, mo, the universal gaugino mass, 
mi/2, the universal soft breaking trilinear coupling constant, Aq, the ratio of Higgs expectation 
values, tan /3, and the sign of the bilinear Higgs coupling constant. In the case of the standard 
cosmology, if we concentrate on smaller values of m and mi/2, then the stau neutralino (fi-%5) 
coannihilation region is the only dark matter allowed region (which is also allowed by the — 2 
constraint) [12J. However, due to the presence of the extra factor R, the parameter space in the 
SSC scenario requires larger values of itiq. This is because a smaller annihilation cross section 
is required in the presence of the dilaton contribution. The magnitude of mo is however, still 
much smaller than the focus point region [13] (in this region, the magnitude of /i is small and 
therefore the lightest neutralino has a large Higgsino component). This difference in mo will 
produce new types of signals at the LHC for the SSC model. For example, in the case of 
the standard cosmology, the allowed region for low mo requires the fi and Xi to have nearly 
degenerate masses within ~ 10 GeV. This produces low energy r's in the final states [HI EE US] . 
In contrast, in the SSC model, Z bosons, Higgs bosons or high energy taus appear in the final 
states. These final states, which we will discuss for this SSC scenario, actually exist in most 
regions of the SUSY parameter space. Therefore, even without any cosmological motivation, 
searching for these signals is a worthwhile exercise. Furthermore, even though we have used 
the SSC as our motivation to probe the signals of the SUGRA model, one can come up with 
any other cosmological framework where the Boltzmann equation is modified in such a way 
that the universe is not really overdosed in this wide region of SUGRA paramater space. This 
analysis is valid for all these scenarios. One can also use other SUSY breaking scenarios in the 
context of SSC. Our analysis of signals can still be applicable in those new scenarios. However, 
additional observables may be required in order to determine the model parameters. 

The determination of the factor R in Eq. [1] is important since it will tell us whether we 
satisfy the cosmological observation for the evolution of dark energy for < z < 1.6. In 
order to calculate R we need to calculate the relic density precisely at the collider. In this 
paper, we first show how to analyze and develop appropriate cuts to extract the signals in the 
newly allowed parameter space in order to determine model parameters. We also construct new 
observables necessary for the determination of such parameters. Then, using these parameters, 
we determine the accuracy of the result for the dark matter content. Finally, when the LHC 
will be operating, the dark matter direct detection experiments also will be probing the SUSY 
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parameter space. The neutralino-proton scattering cross section is different for this newly 
allowed parameter space compared to the standard cosmology case. Thus we also determine 
the accuracy of the result for the neutralino p roton scattering cross section based on the 
LHC measurements. This will be very useful when we will combine the data from these direct 
detection experiments with the LHC data to extract the final model. 

The remainder of this paper is organized as follows: In Section 2, we discuss the parameter 
space of this model and compare with the standard cosmology, followed by characterizing 
the SUSY signals at the LHC in Section 3. In Section 4, we show the determination of model 
parameters and the prediction of relic density and neutralino-proton cross section. We conclude 
in Section 5, where some comments on the applicability of our main results to other string theory 
models are also discussed briefly. 

2 Parameter Space 

The mSUGRA model parameters are already significantly constrained by various experimental 
results. Most important for limiting the parameter space are: (i) the light Higgs mass bound 
of M h o > 114 GeV from LEP [T7|, (ii) the b — > 37 branching ratio bound of 1.8 x 10~ 4 < 
B{B — > Xg'-f) < 4.5 x 10~ 4 (we assume here a relatively broad range, since there are theoretical 
uncertianties in extracting the branching ratio from the data) [18) 1 ! . (iii) the 2a bound on the 
dark matter relic density: 0.095 < QcDMh 2 < 0.1117 |T], (iv) the bound on the lightest chargino 
mass of M~± > 104 GeV from LEP [22] and (v) the muon magnetic moment anomaly a M , where 

the present deviation from the SM value is (29.5 ± 8.1)10- 10 [2D E£2 [21 125] • Assuming the 
future data confirms the a M anomaly, the combined effects of — 2 and > 104 GeV then 

only allows /i > 0. Figure [1] shows g-2 curves for 1.07 x 10 -9 , 1.91 x 10~ 9 , 3.59 x 10~ 9 and 
4.43 x 10—9 which are within two sigma deviation. 

Since the mSUGRA parameters determine the masses of our supersymmetric particles, they 
also determine the dark matter (or neutralino) relic density. We can find which regions of the 
parameter space will agree with the dark matter relic density observed by WMAP. The region 
allowed by WMAP for the standard big bang cosmology is vastly different from the region for 
the SSC model (8]. The comparison of these two regions is shown in Fig.[T]for the case of A = 
and tan (3 = 40. We see a clear separation between the standard big bang cosmology region and 
the SSC region. The standard big bang cosmology region is the very narrow region ranging from 
350 GeV < mi/2 < 900 GeV (850 GeV < M~ g < 2000 GeV), and 200 GeV < m < 350 GeV. 
The SSC region is much broader for m 1( / 2 < 800 GeV, and is higher in mo which ranges from 
400 to 500 GeV. 

1 The present experimental world average for b — » S7 is (3.52 ±0.23 ±9) x 10—4 [19] and the SM contribution 
has been evaluated to be (3.15 ± 0.23) x 10—4 [20] . The b — > S7 constraint does not have much of an impact 
in our analysis. The signals, Higgs+jets+ fir and 2 tau+jets+ ftr in our study are available for a large region 
of parameter space. The Z+jets+ fjj final states however arise in the parameter space where b — > sj is large 
which we mention later. We, however, still discuss this final state since it is easy to evade the b — > sj constraint 
without much change in the final state. 
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Figure 1: WMAP allowed parameter space for the SSC and standard big bang cosmology shown 
for Aq = tan/3 = 40. The very thin green (grey) band is where the neutralino relic density 
calculated by standard big bang cosmology agrees with the WMAP3 limits 0.0950 < Qh 2 < 
0.1117. The thicker dark purple band shows the same agreement using the SSC calculation of 
the relic density. Also shown are the Higgs mass boundary (dashed dotted blue line), muon 
— 2 boundaries (dashed and dotted red lines), a hatched cyan region which is excluded by 
b — > S7 experimental bounds, and a lower solid red region where the neutralino is not the 
lightest supersymmetric particle. 
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Figure 2: The dominant decay branching ratios of decays from the x\ are shown here. Each of 
the four plots shows how the branching ratios vary with m at constant m^. Together they 
survey the SSC band of parameter space in Fig. [TJ 

3 Signals at the LHC 

The SSC region of parameter space has some unique characteristics which are distinguishable 
from the Trx? coannihilation region which appears for the lower values of m . For example, 
let us consider the decay chains of the dominant SUSY production mechanism at the LHC, 
which will produce the squark and gluino, in pairs (e.g., qg). In the coannihilation region, the 
dominant decay chain for the squark (gx) is qi — > qx® ~^ QT^i Q^Xi- Here X2 * s the second 
lightest neutralino. Thus, this region produces r's, along with jets and missing transverse 
energy, fa. However, the characteristic decay in the SSC region is X2 ~^ ^°Xi- I n this case, 
we would expect h° — > bb, along with jets and fa. Figure [2] shows the branching ratios for the 
X2 decay. As we increase rai/2, the branching ratios shift from Higgs dominant decay chains to 
t dominant decay chains. However, we will easily distinguish this SSC r dominated region fr 
om the coannihilation region by observing a large mass difference between the f\ and x?- For 
even lower values (m^ < 350 GeV), the X2 ~^ ^Xi decay becomes dominant. Typical 
mass spectra are shown for points in the Higgs boson, Z boson and two r final state dominated 
regions in Tables [Q, [2j and El 

3.1 Signals 

We have three possible signals in this model. These signals are the following: 

• Higgs + jets + fa 

• Z + jets + fa 



5 



Table 1: SUSY masses (in GeV) and dominant branching ratios for x 2 f° r the point itlq = 

471 GeV, rui/j = 440 GeV, tan (3 = 40, A = 0, and fi > 0. Notice that the x 2 ~> f i r is 

kinematically forbidden. For this point, fl~oh 2 = 0.089 and o ~o = 2.42 x 10~ 9 pb. The total 

Xi p-Xi 
production cross section for this point is a — 1.61 pb. 
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Table 2: SUSY masses (in GeV) and dominant branching ratios for x 2 f° r the point itlq = 
471 GeV, mi/2 = 320 GeV, tan/3 = 40, A = 0, and fi > 0. We chose this point to examine 
despite the fact that it is within the region excluded by b — > S7. We did this to examine the 
behavior of X2 ~^ ^Xi & t its maximal branching ratio. The total production cross section for 
this point is a = 7.10 pb. 
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Table 3: SUSY masses (in GeV) and dominant branching ratios for x 2 f° r the point mo = 

440 GeV, m V 2 = 600 GeV, tan/3 = 40, A = 0, and \i > 0. For this point, fi~ /i 2 = 0.106 and 

Xi 

0" _^o = 7.19 x 10~ 10 pb. The total production cross section for this point is a = 0.446 pb. 
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• 2 r + jets + $i 
3.1.1 Higgs + Jets + tfk 

The Higgs + jets + |r signal appears in the lower myi region (400 < my2 ^ 500 GeV) within 
the SSC band of parameter space in Fig. [TJ The Higgs + jets + |r signal is characterized by 
the decay chain, qi — > qh°X\- The Xi does not interact in the detector, and thus leaves 

a large |r signal. The h° and jet carry information about the SUSY particles in this chain. In 
particular, the h° + jet invariant mass distribution has an endpoint which depends upon their 
masses: 




(2) 

For our analysis, we generate events using PYTHIA [27], which is linked with ISASUGRA [28] to 
generate the mSUGRA mass spectrum. We pass these events to a detector simulator called 
PGS4 [29]. 

To measure the endpoint, we begin by selecting our events with the following cuts |15j : 

• At least two jets with px > 200 GeV as well as \rj\ < 2.5, 

• |r > 180 GeV, 

• Pt 1 + Pt* 2 + It > 600 GeV, and 

• at least two 6-tagged jets [29] with pr > 100 GeV and \r}\ < 1.5. 

These cuts remove the majority of the SM background, such as tt, W+jets, and Z+jets [26], as 
well as some background from SUSY events which do not contain the decay chain q^ — > qx\ — * 
qh Xi- Unforeseen SM backgrounds at the LHC can be shape analyzed from the data and then 
subtracted from our desired signal. 

Next we identify Higgs bosons in the event. We select all pairs of 6-tagged jets with p\ > 
100 GeV and 0.4 < AR bb < 1. The lower AR limit is due to a jet clustering cone size 
in PGS4, while the upper limit is motivated by a study of Higgs decays in mSUGRA events 
(mo = 471 GeV, = 400 GeV) at the generator level. (See Fig. [21) We then form the b 
pair invariant mass. Figure H] shows a peak between 100 GeV and 120 GeV, consistent with the 
Higgs mass, along with a continuum background. For each b pair, we form 2b + jet systems, 
using the two jets with the greatest transverse momenta of the event. These two leading jets 
will primarily come from three different decay chains: 

• q -»• qg 

• qn^ qx°i 

• qL^ qxl -> q h °x°i 
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Figure 3: Correlation between p™ m = mm(p\,, pj) and Ai?^ from h° — > bb at the generator level 
for mSUGRA events at m = 471 GeV, m\/2 = 400 GeV. Here is a separation between 

b and b in 7/-0 space. The inset histogram is the distribution for p™ m > 100 GeV. This 

shows that bb pairs from a single Higgs decay will most often have a separation of AR < 1 for 
6-jets with transverse momentum greater than 100 GeV. Any b pairs not from a single Higgs 
decay will instead have no particular separation. 



Each b pair combined with these two leading jets will form two effective masses, and 
Mbbj 2 . If we combine the jet from qn — > qxi with our Higgs, it can have a larger Mbbj than the 
endpoint expected from the q~i — * QX2 ~* <?^°Xi decay chain. Thus, we simply choose the lesser 
of Mbbj! and M^ 2 , denoting it as M^J d . This selection is a similar to that shown in Ref. [26] . 

At this stage, we still suffer from bb combinatoric background as seen in Fig. HI To estimate 
the combinatoric background in the Higgs mass window, we perform a sideband subtraction 
method. We form the M^ d distribution using b pairs in the Higgs peak window and using 
b pairs from two sideband windows (70-90 GeV and 130-150 GeV) in Fig. HI This second 
distribution is scaled by the ratio of the background shape evaluated in the Higgs window to 
the sideband windows. (See Fig. @J. Then we subtract this scaled sideband M^ d distribution 
from the Higgs M^ d distribution. Since the kinematical endpoint occurs when the Higgs is 
back to back with the jet, we select events with AR^oj > 1.2. 

In order to determine the endpoint, we fit the mass distribution to a combination of a 
Landau probability distribution function, P L , and a straight line: 

f/\ J kP]j\X) %>peaki O") if X <C Xp ea k /cy. 

J\ x ! — \ , p / \ , / \ - f ' \ 6 ) 

I rCiLyXyXpeakytJ) -p CH\X Xp ea fcj II X S> Xp ea k 

where x corresponds to the Higgs plus jet invariant mass, x pea k is the most probable value of the 
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M bb (GeV) 

Figure 4: The invariant mass distribution of PGS b jet pairs. The central blue (dark grey) bins 
are the Higgs peak window. We perform a background subtraction by selecting the green (light 
grey) sideband windows. The background fit is the black curve. 

Landau distribution, k scales the height of the function, and a is the slope of the linear portion. 
Figure [5] shows the fittings at two mSUGRA points around our reference point described in 
Table [IJ One can see the change in the shape and the end point as increases. The 

slight change in shape between the two histograms in Fig. [5] is due to the fact that the SUSY 
background for this signal does not shift as we vary m^, and that it dies off around 800 GeV. 
Thus, the mi/2 = 400 GeV histogram which has an endpoint around 750 GeV has a slight 
shoulder after the endpoint, whereas the mi/2 = 480 GeV histogram has an endpoint around 
900 GeV with no shoulder. 

Despite such shoulders, the shape of the distribution stays similar as we vary the model 
parameter m 1( / 2 . Thus we can use the same fitting function for such points. Also, since the 
endpoint, does not depend on any third generation sparticles (see Eq. [2]), it is independent 
of the parameters A and tan (3. However, if we increase m the situation changes. For higher 
m values the q becomes significantly heavier than the g. The result of this, for instance in the 
case of m = 651 GeV and my 2 = 440 GeV, is that B(q L — > qg) = 10% and B(q R — > qg) = 22%. 
The quark jets from such decay chains are much softer than background jets from lower mo 
points. Thus, the nature of the background changes, which changes the shape of the Higgs plus 
jet invariant mass distribution: The endpoint becomes very sharp. Thus a simple linear fit is 
sufficient to find the endpoint. A sample fit of this higher mo region is shown in Fig. [6j 
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Figure 5: The Higgs (tagged b jet pair) plus jet invariant mass distribution reconstructed 
through PGS in two 500 fb" 1 mSUGRA samples at (m , m 1/2 ) = (471 GeV, 400 GeV) and 
(471 GeV, 480 GeV) for the black histogram with red (gray) fit and blue (gray, filled) histogram 
with dark blue(dark gray) fit, respectively We fix tan/? = 40, A = 0, and fi > 0. 




200 400 600 800 1000 1200 1400 1600 1800 2000 

M 2 b n b f (GeV) 

Figure 6: The Higgs (tagged b jet pair) plus jet invariant mass distribution reconstructed 
through PGS in a 500 fb" 1 mSUGRA sample at m = 651 GeV, m 1/2 = 440 GeV, tan/3 = 40, 
A = 0, and ji > 0. 
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Figure 7: The Z plus jet invariant mass distribution reconstructed through PGS in one 50 fb 1 
mSUGRA sample at m = 471 GeV, m 1/2 = 320 GeV, tan/3 = 40, A = 0, and ji > 0. 



3.1.2 Z + Jets + fa 

The final state of Z + jets + j£j events becomes a key signal in a lower mi/2 region (muj — 
300 GeV) where the x° ~~ * <?^°Xi decays are kinematically suppressed. The decay chain and 
endpoint equation (Eq. [2]) are exactly the same under the replacement of the Higgs boson mass 
with the Z boson mass. 

To construct the Z plus jet invariant mass and measure the endpoint, we follow the very 
same procedure as the Higgs plus jet analysis, but with Z — > 11 decays. We select events with 
the same initial cuts as in the Higgs plus jet analysis and reconstruct the Z — > 11 decays instead 
of the Higgs decays. To select our Z bosons we find pairs of isolated leptons with p T > 20 GeV 
in our event. We keep pairs of leptons with invariant mass within the Z mass window, where 
85 GeV < Mu < 97 GeV. Then we form 21 + jet systems using the two jets with the greatest 
transverse momenta. We again keep only the lesser of the two 21 + jet invariant masses, Mf^ d . 
To ensure we select mostly Z bosons within this signal we use an opposite-sign-same-flavor 
minus opposite-sign-opposite-flavor subtraction. A sample distribution of the result is shown 
in Fig. [71 This figure shows us a measurable endpoint very similar to that of the Higgs plus jet 
invariant mass technique. 

3.1.3 2r + Jets + $r 

The 2r + jets + fa signal appears in the higher myz region (mi/ 2 > 500 GeV) within the SSC 
band of parameter space in Fig. [TJ The full decay chain qi — > 9X2 ~^ ~^ rT Xi produces a 
characteristic final state consisting of r's, high pt jets from the q~L decay, and fa from the Xv 
Again, the two r particles and the jet carry information about the supersymmetric particles 
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in the decay chain. The visible ditau invariant mass distribution and the 2r plus jet invariant 
mass distribution both have endpoints depending on the supersymmetric particle masses: 
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We again make use of such kinematical observables by measuring these endpoints. However, in 
this case, we are restricted by the background for M^ d . To avoid this background, we measure 
the peak instead, since the peak is proportional to the endpoint. Events were generated using 
ISA JET [28J and the detector effects were simulated using PGS4 [29|. 

To measure these observables, we select our events with the following cuts [14J: 

• At least two jets with > 200 GeV as well as \rj\ < 2.5, 

• fir > 180 GeV, and 

• p^t 1 + P>t 2 + fir > 600 GeV. 

• We reject events where either one of the two leading jets is tagged as a b jet [29|. 

These cuts remove the majority of SM backgrounds (tt, W+jets, and Z+jets), as well as 
background from SUSY events containing stops or sbottoms. Here again, unforeseen SM back- 
grounds at the LHC can be shape analyzed from the data and then subtracted from our desired 
signal. We do not discuss the details of the event selections, but instead refer the reader to our 
previous and ongoing studies [15], [161 E] • 

Sample M TT and Mj TT distributions for points similar to that shown in Table E] are displayed 
in Figs. [H] and O These figures also show how the peak and endpoint shift under changes of 
mi/2. Since the endpoint of the 2r + jet invariant mass distribution, Mj^, does not depend on 
any third generation superparticles (see Eq. EJ), it will only shift for variations of the mSUGRA 
parameters m and mi/ 2 . However, the peak of the 2r invariant mass distribution, M^ ak , 
depends on the stau mass (see Eq. HJ), and will thus depend on all four mSUGRA parameters. 



4 Determining Model Parameters 

We have shown in previous works [T5J EE] that we can obtain mass measurements of the su- 
persymmetric particles in the neutralino-stau coannihilation region by utilizing each final state 
and parameterizing kinematical observables, such as those described in the previous section, in 
terms of the SUSY masses. Our goal is to determine the mSUGRA model parameters mo, m\/2, 
Aq, and tan/3 since we want to determine the dark matter content and the neutralino-proton 
cross section. The fifth mSUGRA model parameter, sign(/i), is assumed to be positive, since 
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Figure 8: The ditau invariant mass distribution reconstructed through PGS in two 500 fb _1 
mSUGRA samples at (m , m 1/2 ) = (440 GeV, 625 GeV) and (440 GeV, 575 GeV) for the 
black and blue(gray, filled) histograms, respectively. We fix tan/? = 40, A = 0, and /i > 0. 

this is preferred by measurements of the b — > 57 branching ratio and the muon g M — 2 [23J. 
To determine the mSUGRA parameters, we will thus need four kinematical observables which 
are linearly independent functions of those parameters. The determination of the parameters 
is then accomplished by inverting four such functional relationships. 

As discussed above, certain regions of the mSUGRA parameter space might give rise to very 
different signals. For each region we can combine different observables to determine the four 
mSUGRA parameters. However, so far, there are not four observables which can be made for 
each signal described above. Thus we introduce the following additional kinematical observables 
which are valid in any region: M e g, , and Mjy . 

The effective mass, M e g, is defined by 

M eff = ^ + p$ t2 + # 3 + $ u + fa, (6) 

where all four leading jets are not 6-tagged jets. This combination carries the information of the 
characteristic SUSY scale. The majority of the of the jets is characteristic of the gluino and 
first two generation squark decays, and the majority of the fa is due to the lightest neutralino 
escaping the detector. As such, the observable M e g depends only on the mSUGRA parameters 
mo and mu2- This is because the parameters Aq and tan/3 only affect the third generation 
superparticles. 

When we construct the effective mass distribution, we use the following cuts[26j: 

• At least one jet with pr > 100 GeV and an additional three jets pt > 50 GeV, where all 
such jets have \r)\ < 2.5, 
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Figure 9: The 2r plus jet invariant mass distribution reconstructed through PGS in two 500 fb 
mSUGRA samples at (m , m 1/2 ) = (440 GeV, 625 GeV) and (440 GeV, 575 GeV) for the 
black histogram with red(gray) fit and blue(gray, filled) histogram with dark blue(dark gray) 
fit, respectively. We fix tan/3 = 40, A = 0, and /i > 0. 



• No isolated leptons with \r]\ < 2.5, 

• fir > 200 GeV and fir > 0.2 x M efr , 

• Transverse Sphericity, St > 0.2, 

• We reject events where any of the four leading jets is tagged as a b jet [29J . 

We find the peak value with an iterative fitting technique. First, we fit the distribution itera- 
tively with an asymmetric gaussian function. The purpose of the iterative fit is simply to find 
the ideal fitting range. Once that is found, we fit with a cubic polynomial to find the peak 
position. A sample effective mass distribution showing the result of this procedure is shown in 
Fig. [[0j This figure also shows that as mi/ 2 increases, the peak increases. 

Two very similar observables can also be constructed. The b effective mass, M^, is defined 

by 

and the 2b effective mass, M^g , is similarly defined by 

M e ( f = p5 tl(6) + + # 3 + ^ + fir. (8) 

Here there are no restrictions on the non-leading jets; they can be either 6-tagged or not. By 
including the leading b jets, which are primarily decay products of the superpartners to the 
third generation quarks, we include information about the parameters A Q and tan f3. 
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Figure 10: The effective mass distribution reconstructed through PGS in two 500 fb _1 mSUGRA 
samples at (m , m l/2 ) = (471 GeV, 400 GeV) and (471 GeV, 480 GeV) for the black histogram 
with red(gray) fit and blue(gray, filled) histogram with dark blue(dark gray) fit, respectively. 
We fix tan (3 = 40, A = 0, and fi > 0. 



To construct these distributions, we use the very same cuts as we used for M eff , with some 
exceptions. For the distribution, the leading p T jet must be tagged as a b jet, otherwise 
we reject the event. For the distribution, the two leading jets must both be tagged 

as b jets. Again, the non-leading jets in and can be either 6-tagged or not. We also 

use the same fitting algorithm for these distributions as we have used for M e g. Sample 
and M e g 6 ' ) effective mass distributions are shown in Figs. [TT] and [T^J 

Now that we have all the observables we need, we can determine the four mSUGRA param- 
eters in any region. We will now describe examples of this method for a Higgs boson dominated 
region, a f dominated region, and a Z boson dominated region. 

Higgs + jets + f)j<: For the Higgs dominant region, we use the following four observables 
to determine our mSUGRA parameters: 

• Effective Mass: = fi(m , m ly / 2 ) 

• b Effective Mass: M^q peak = f 2 (rn , rni/ 2 , A , tan f3) 

• 2b Effective Mass: M^ 1 peak = f3(m , mi/ 2 , A , tan/3) 

• Higgs plus jet invariant mass: M^ d,end = /4(m ,mi/ 2 ) 

These functional forms are determined by examining how each kinematical observable changes 
while varying one of the mSUGRA parameters. Examples of this are shown in Figs. [T3land[T4l 
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Figure 11: Same as Fig. [TDJ except that this is the b effective mass distribution and that the 
peak fits are not shown. 



To determine our mSUGRA parameters, we invert these functional forms into functions of the 
mSUGRA parameters in terms of the kinematical observables. Then we can simply plug in the 
values of the observables into the inverted functions to solve for the mSUGRA parameters. To 
get the uncertainties of the mSUGRA determinations, we propagate the uncertainties of the 
measured observables through the inverted functions using a Monte Carlo method. 

We perform a sample analysis for the Higgs region with the following result: mo = 472 ± 
50 GeV, rrii/2 = 440 ±15 GeV, A = 0±95 GeV, and tan/? = 39 ±17. These uncertainties were 
achieved at 1000 fb~ . The relation between the uncertainties and the luminosity is shown for 
these parameters in Figs. [15] and [16j Using these results, we can also calculate the neutralino 
relic density and proton-neutralino cross section. The result for 1000 fb" 1 is: tt~ h 2 = 0.10 ± 

0.15 and a ~o = (1.9±3.7) x 10~ 9 pb. The uncertainty ellipses on the f2~o/i 2 -tan/3 and a ~o- 
p-Xi Xi p-Xi 

tan (3 planes are shown in Fig. [171 Since the uncertainties in each of these values are larger than 

100%, the uncertainty ellipses get pushed into negative (unphysical) values of Q~oh 2 and a ~o. 

Xi P~ Xi 

As such, we have cut these ellipses off at the x-axes in Fig. effigError Ellipse: HiggsRegion. 

Z + jets + ft?: The analysis technique in the Z dominant region is just the same as 

the Higgs region if we replace the Higgs plus jet invariant mass with the Z plus jet invariant 

mass. The endpoint of the latter can be measured with better precision. This is due to both an 

increase in production cross section and the ease of reconstructing Z bosons from lepton pairs. 

This results in a more precise determination of Q~oh 2 and a ~o comparable to that of the f 

Xi p-Xi 

dominant region shown below. However, we suffer from small B (Z — > 11) values. For a useful 
measurement, we need E (x° ~~ > ^Xi) ~ 50%. However, such a region does not exist outside 
of the b — ► S7 bound, as shown in Fig. [181 For lower tan f3, the same conclusion holds since 
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Figure 12: Same as Fig. [TT] except that this is the 2b effective mass distribution. 



we get constraints on the smaller values of mi/ 2 due to Higgs mass. Therefore, we do not go 
into detailed analysis of the determinat ion of model parameters in the Z + jets + fir region. 
However, one can use the observables of the Higgs + jets + fir region, e.g., M e g, , and 
M^ b \ to reconstruct the model parameters. 

2 t + jets + fir: For the f dominant region, we use the following four observables to 
determine our mSUGRA parameters: 

• Effective Mass: Mjg ak = /i(m ,mi/ 2 ) 

• b Effective Mass: M^q peak = f 2 (m 0} mi/ 2} A ,tan f3) 

• Ditau Invariant Mass: M^ ak = f 3 (m , m 1 / 2 , A , tan j3) 

• Ditau plus jet invariant mass: Mj^ k = ^(mo,™^) 

We again perform an inversion to determine our mSUGRA parameters, as well as propagate 

the uncertainties in the same way as in the Higgs region. Our sample analysis for this region 

yields: m = 440 ± 23 GeV, m 1/2 = 599.9 ± 6.1 GeV, A = ± 45 GeV, tan/? = 40.0 ± 2.7, 

Q~oh 2 = 0.103 ± 0.019, and a ~o = (7.6 ± 1.6) x 10~ 10 pb. These uncertainties were achieved 
Xi p-Xi 

at 500 fb~ . The parameter tan/? is determined with much higher accuracy since we can use 

observables involving staus; the staus are very sensitive to tan/3. Again we show the relation 

between the uncertainties and the luminosity for this result in Figs. [19] and [20l We also again 

show the uncertainty ellipses on the Q ~o/i 2 -tan 8 and o ~o-tan 3 planes in Fig. EU Since tan 3 

Xi p-Xi 

is determined with better accuracy compared to the Higgs dominant region, the relic d ensity 
and proton-neutralino cross section are also determined with a better accuracy. 
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Figure 13: The left plot shows the change in M^ d ' end under variations of m (mi/2 = 440 GeV) 
within the Higgs dominant region of parameter space. The right plot shows the same for 
variations in mi/2 (mo — 471 GeV). Combining the functions plotted results in the functional 
form M^ d ' end = / 4 (m ,m 1; /2). The la uncertainty bands (dashed lines) represent 500 fb _1 of 
data. 
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Figure 14: The left plot shows the change in M eff pea under variations of Aq within the Higgs 
dominant region of parameter space (mo = 471 GeV, mi/2 = 440 GeV). The right plot shows 
the same for variations in tan (3. Combining the functions plotted along with variations in mo 
and mi/2 results in the functional form M e g peak = ^(mo, m^, Aq, tan (3). The la uncertainty 
bands (dashed lines) represent 1000 fb _1 of data. 
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Figure 15: The left plot shows the change in the measurement uncertainty of m in the Higgs 
dominant region of parameter space (m = 471 GeV, mi/ 2 = 440 GeV) for different luminosities. 
The right plot shows the same for the uncertainty in m 1( / 2 . 
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Figure 16: Same as Fig. [15] but for A Q and tan/3. 
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Figure 17: The left plot shows the la uncertainty ellipse on the Q~oh 2 -ta,nf3 plane in the Higgs 

Xi 

dominant region of parameter space (mo = 471 GeV, m-1/2 = 440 GeV). The right plot shows 

the same for the a ~o-tan/3 plane. The uncertainty in tan/3 is the main source of both the 
p-Xi 

uncertainty of Q~oh 2 as well as the uncertainty of a ~n. These results are for 1000 fb _1 of 

Xi p-Xi 

data. 
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Figure 18: The branching ratios for X2 ~^ h Xi an d X2 ~^ ^Xi as a function of mi/2 for 
m = 471 GeV, A Q = 0, and tan f3 = 40. Also shown is the b — > sj exclusion region (cyan filled 
region) from FigJTJ 
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Figure 19: Same as Fig. [15] except within the r dominant region of parameter space (mo 
440 GeV, m 1/2 = 600 GeV). 
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Figure 20: Same as Fig. [15] except for A and tan/? within the r dominant region of parameter 
space (m = 440 GeV, m 1/2 = 600 GeV). 
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Figure 21: Same as Fig. [T7] except within the f dominant region of parameter space (mo : 
440 GeV, m 1/2 = 600 GeV). These results are for 500 fb" 1 of data. 
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5 Conclusions and Discussion 



In this paper, we have studied the mSUGRA final states at the LHC which are motivated 
by supercritical string cosmology (SSC). In the SSC case, the time dependent dilaton not 
only contributes to the dark energy but also to the Boltzman equation which determines the 
dark matter content of the universe. Consequently the dark matter profile in this model is 
different compared to the standard cosmology. We found that the dark matter allowed region 
has larger values of m compared to the standard cosmology case. Thus, the final states in the 
SSC scenario are different from those of the standard cosmology. For example, in the case of 
standard cosmology for smaller values of m (also allowed by the — 2 constraint), we have 
low energy taus in the final state due to the proximity of the stau to the neutralino mass in 
the stau-neutralino coannihilation region. On the other hand, in the SSC case the final states 
contain Z bosons, Higgs bosons or high energy taus. In fact these final states dominate in most 
of the allowed SUGRA parameter space. Therefore, by analyzing the parameter space of the 
SSC model we actually investigate most regions of the SUGRA parameter space at the LHC. 
We analyzed the signals involving Higgs + jets + fir, Z + jets + fir, and 2r + jets + fir and 
constructed observables such as the endpoints of invariant mass distributions M^, Mzj, and 
Mj TT and the peak point of M TT . 

In order to determine all parameters of the mSUGRA model we needed additional obsevables 
such as the peak positions of the effective masses M^f and M^ b) . These observables are used 
for determining tan/3. We found that m , mi/2, and tan/5 can be determined with 11%, 3% 
and 44% accuracies respectively in the Higgs boson dominated final states region for 1000 fb _1 
of data. The Z boson dominated final state region is mostly ruled out by other experimental 
data. However, the technique used to analyze the Z boson dominated region is nearly identical 
to that of the Higgs boson dominated region. In the stau dominated region, m , mi/ 2 , and 
tan/5 can be determined with 5%, 1% and 7% accuracies respectively for 500 fb _1 of data. 
The accuracy of determining tan/3 is improved in the tau dominated final state region since 
we use observables involving the staus which are very sensitive to the variation of tan/3. Once 
all the parameters are known, the dark matter content can be determined in all these cases. 
In the Higgs dominant case, the accuracy of determining the dark matter content is 150% for 
1000 fb -1 of data. In contrast, the accuracy of relic density in the stau dominated region is 
18% for 500 fb _1 of data, which is much better due to a higher accuracy of tan (3 determination. 
These techniques can be applied in the case of nonuniversal models as well, where we will need 
more observables to determine the model parameters. 

When the LHC will be operating, we will also have results from the dark matter direct 
detection experiment. We found that the cross section for these experiments can be predicted 
from the LHC measurements with an accuracy of 195% for 1000 fb" 1 of data in the Higgs boson 
dominated region and 21% for 500 fb~ x of data in the stau dominated region. This cross section 
however includes uncertainity due to the form factors. 

As a remark, our phenomenological study assumed several key detector performances of the 
present ATLAS and CMS detectors, such as b-tagging and tau identification efficiencies. We 
find that one needs 500-1000 fb" 1 of data. The regime of such high luminosity can be realized 
with the LHCs luminosity upgrade as well as the upgrade of both ALTAS and CMS detectors. 
Thus, our results are just a guideline for the physics case if the performance of both upgraded 
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detectors is the same even at such high luminosity operation of the luminosity-upgraded LHC. 

In this analysis we examined the overdense region of the mSUGRA model since the underly- 
ing cosmological theory converts the overdense region into a region with correct relic abundance. 
This analysis holds for any cosmological model with similar features. 

Before closing we repeat some cautionary remarks regarding the microscopic model depen- 
dence of such studies [11] . As already mentioned in the introduction, the low-energy limit of 
string theory is incredibly non-unique, as it depends on the complicated details of compactifi- 
cation and SUSY breaking procedures. Various models lead to different predictions, and some 
of them may lead to completely different phenomenology as far as dark matter studies are con- 
cerned. For instance, there are heterotic string models entailing non-thermal dark matter [10J, 
whose detection requires totally different techniques from the ones employed here. 

Nevertheless, there are string models which can be analyzed rather generically within the 
methods outlined in this work, in the sense that the observables discussed in this analysis can 
also be used to extract information on dark matter in such string-inspired models as well. For 
instance, the moduli-dominated sector of the heterotic (orbifold-compactified) class of models 
examined in [10] has five parameters the gravitino mass 7B3/2, the vacuum expectation value of 
the real part of the (uniform) Kahler modulus (t + t), the modular weights of the Pauli- Villars 
regulators parameterized by p, the value of the Green-Schwarz coefficient Sqs an d tan/3. The 
parameters of this model can be determined in the same spirit as shown in the paper and 
thereby the dark matter density can also be determined in the way we have described in 
this work. The same procedure can be applied to the dilaton dominated models described in 
the same reference. Depending on the model which is used, we may need to construct more 
observables to determine all the model parameters. 
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